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Abstract 

Background: The alternative quantum mechanical description of total energy given by Bohmian theory was 
merged with the concept of the golden ratio and its appearance as the Heisenberg imbalance to provide a new 
density-based description of the valence atomic state and reactivity charge with the aim of clarifying their features 
with respect to the so-called DFT ground state and critical charge, respectively. 

Results: The results, based on the so-called double variational algorithm for chemical spaces of reactivity, are 
fundamental and, among other issues regarding chemical bonding, solve the existing paradox of using a cubic 
parabola to describe a quadratic charge dependency. 

Conclusions: Overall, the paper provides a qualitative-quantitative explanation of chemical reactivity based on 
more than half of an electronic pair in bonding, and provide new, more realistic values for the so-called "universal" 
electronegativity and chemical hardness of atomic systems engaged in reactivity (analogous to the 
atoms-in-molecules framework). 

Keywords: Electronegativity, Chemical hardness, Bohmian mechanics, Heisenberg imbalance equation, Slater 
electronic density 



Introduction 

Recently, the crucial problem regarding whether chem- 
ical phenomena are reducible to physical ones has had an 
increasingly strong impact on the current course of con- 
ceptual and theoretical chemistry. For instance, the fact 
that elements arrange themselves in atomic number (Z) 
triads in approximately 50% of the periodic system seems 
to escape custom ordering quantifications [1,2]. The same 
applies to the following: the fascinating golden ratio (r) 
limit for the periodicity of nuclei beyond any physical 
first-principle constants, which provides specific periodic 
laws for the chemical realm [3-6]; the fact that atoms 
have no definite atomic radii in the sense of a quantum 
operator, and even the Aufbau principle, which, although 
chemically workable, seems to violate the Pauli Exclusion 
Principle [7]; at the molecular level, the well-celebrated 
reaction coordinate, which, although formally defined in 
the projective energy space, does not constitute a variable 
to drive optimization in the course of chemical reactions, 



Correspondence: mv_putz@yahoo.com 

Laboratory of Computational and Structural Physical Chemistry, Biology- 
Chemistry Department, West University of Timisoara, Pestalozzi Street No. 16, 
Timisoara RO-300115, Romania 



appearing merely as a consequence of such reactions [8]; 
the problem of atoms in molecules [9], i.e., how much of 
the free atoms enter molecules and how much indepen- 
dency the atoms preserve in bonding; and chemical 
bonding itself, which ultimately appears to be reinter- 
preted as a special case of bosonic condensation with the 
aid of bondons - the quantum bosons of chemical bond- 
ing, which, without being elementary, imbue chemical 
compounds with a specific reality [10,11]. 

In the same context, the specific measure of chemical 
reactivity, electronegativity (x), which lacks a definite 
quantum operator but retains an observable character 
through its formal identity with the macroscopic chem- 
ical potential x = ~H [12,13], was tasked with carrying 
quantum information within the entanglement environ- 
ment of Bohmian mechanics [14-17] and has thus far 
been identified with the square root of the so-called 
quantum potential % = Vq 2 [6] . 

However, the striking difference between an atom as a 
physical entity, with an equal number of electrons and 
protons (thus in equilibrium), and the same atom as a 
chemical object, with incomplete occupancy in its per- 
iphery quantum shells (thus attaining equilibrium by 
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changing accepting or releasing electrons), is closely 
related to the electronegativity phenomenology in model- 
ing chemical reactivity. Moreover, this difference triggers 
perhaps the most important debate in conceptual chem- 
istry: the ground vs. valence state definition of an atom. 

The difficulty may be immediately revealed by consider- 
ing the variation in the total energy (of the ground and/or 
valence state - see below for an explanation of their dif- 
ference) around the physical equilibrium (neutral atom) 
attained between the release (by ionization, I) and receipt 
(through affinity, A) of electrons toward chemical equi- 
librium (in molecules, chemical bonding). Accordingly, 
the curve passing through these points apparently only 
behaves as shown in Figure 1(a), while in all systems 
(with numerical / and A), the obtained interpolating curve 
presents a minimum toward accepting electrons (see 
Figure 1(b)), thus confirming the electronegativity concept 
as a chemical reality, although with a predicted fractional 
charge (for example, the critical charge N ) on an atom at 
chemical equilibrium (i.e., not reducible/comprehensible 
to/by an ordinary physical description of atoms). 

However, the physical-to-chemical paradox continues 
in an even more exciting fashion as follows. When, in 
light of the above discussion, electronegativity is recog- 
nized with the two-point limits shown in Figure 1(b), 
namely [13,18] 

[I ...AW e [-1,0) 

\A...N react e (0,+l] 1 ; 

the limits represent tangents to a curve that does not 
describe chemical equilibrium but an excited state driven 
by the parabolic form 

Edft = -xN + IN 2 (2) 



which happens to correspond to the celebrated den- 
sity functional theory (DFT) working energy expression 
[13,19-21] written in terms of electronegativity and chem- 
ical hardness, respectively defined as follows [13,22-25]: 

=4 (II) vw = Kli),^ (4) 

The point is that curve (2) is not chemically minimized, 
although it is very often assumed to be in the DFT 
invoked by the chemical reactivity literature [13,26-29]; 
however, the curve cannot be considered indicative of a 
sort of ground state (neither reactive nor critical states of 
Figure 1). Additionally, by comparing the curves of 
Figure 1 (a) and (b), the curve of eq. (2) occurs above both 
the reactive and critical curves of Figure 1; it thus should 
represent the chemical valence state with which to oper- 
ate. Therefore, much caution should be taken when work- 
ing with eq. (2) in assessing the properties of atoms, 
molecules, atoms in molecules, etc. Nevertheless, this is 
another case of chemistry not being reducible to physics 
and should be treated accordingly. It is worth noting that 
Parr, the "father" of eq. (2) and a true pioneer of concep- 
tual density functional theory [30,31], had tried to solve 
this dichotomy by talcing the "valence as the ground state 
of an atom in a perturbed environment". This statement is 
not entirely valid because perturbation is not variation 
such that it may be corrected by applying the variational 
principle to eq. (2), for example. In fact, using such 
variation should be considered a double variational tech- 
nique that is necessary to arrive at the celebrated 
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chemical reactivity principles of electronegativity and 
chemical hardness, as recently shown [32]. 

The current line of work takes a step forward by 
employing the double variation of the parabolic energy 
curve of type (2) to provide the quantum (DFT) valence 
charge of an atom (say, N ) and to compare it either 
quantitatively and qualitatively with the chemical critical 
charge N . The goal of these efforts is to gain new 
insight into the valence state and chemical reactivity at 
the quantum level. To this end, the relation of Bohmian 
mechanics to the concept of the golden ratio will be 
essential and will be introduced in the following. 

The consequences of the joint consideration of Bohmian 
mechanics and the golden ratio for the main atomic sys- 
tems will be explored, and the quantum chemical valence 
state will be accordingly described alongside the so-called 
universal electronegativity and chemical hardness, refining 
the work of Parr and Bartolotti [33] as well as generalizing 
the previous Bohmian-Boeyens approach [3,4]. 

Background methods 

Two apparently disjoint theories of matter will be 
employed to characterize the quantum valence of an 
atom: Bohmian mechanics - furnishing the main equa- 
tion for total energy - and the fundamental quantum 
mechanics through the Heisenberg combined with de 
Broglie principles providing the wave-particle indeter- 
minacy framework in which the golden ration depend- 
ency of Z/N naturally appears as quantifying the valence 
states of atoms considered the "ground state" of the 
atomic chemical reactivity. 

Bohmian mechanics 

Because of the need to reduce Copenhagen's indeter- 
minacy for quantum phenomena, i.e., by associating it 
the quantum description of "Newtonian" forms of motion, 
though by preserving probability densities, quantum 
averages, etc., the so-called "minimalist" quantum theory 
may be formulated following the Bohm quantum mech- 
anical program as follows. 

One begins with the general eikonal wave-function 
form [14] 



f{r,t)=R{r) exp[^S{r.,t) 



(5) 



which represents the mid-way between wave and particle 
mechanics because it contains both information regard- 
ing Hamilton-Jacobi theory and the Wentzel-Kramers- 
Brillouin (WKB) approximation [34] through the principal 
phase function S{r,t) while preserving the amplitude 
relationship with the systems' quantum density: 



In this framework, the Schrodinger equation, 

3 h 2 
ih - f(r, t) = ~ — V 2 Mr, t) + V(r)y/(r, t) (7) 

decomposes into real and imaginary parts. The real part 
can be expressed as follows: 



dS(r,t) (V r S(r,t)) 2 h 2 V 2 r R{r 



dt 



2m 2m R(r) 



+ V(r) = 0 (8) 



representing a continuous "fluid" of particles driven by 
the "guidance" momentum: 



mv = p = V r 5(r, t) 



(9) 



moving under a joint external potential V(r) as well as 
under the so-called quantum potential influence: 



fV V 2 R{r) 
2m R(r) 



(10) 



The consequences are nevertheless huge. For example, 
this methodology allows for the interpretation of the tra- 
jectories orthogonal to constant surfaces, by cancelling 
the Laplacian of the wave fronts V r S(r, t) = 0, which are 
obtained from eqs. (8) and (9) as the quantum equation 
of motion: 



dp 
dt 



-V r [V Q (r) + V(r)\ 



(11) 



Equation (11) resembles the classical Newtonian 
acceleration-force relationship only in a formal way; in 
fact, it generalizes it: it prescribes acceleration motion 
even in the absence of an external classical potential. 
This is essential in explaining why the inter-quark forces 
increase with the increase in inter-quark distances, no 
matter how great a separation is considered (a specific 
quantum effect), due to the presence of a quantum po- 
tential that does not fall off with distance as V does. It 
also nicely explains the observed interference patterns in 
double-slit experiments in the absence of classical forces. 
Alike, eq. (11) also appears suited for modeling chemical 
reactivity for the valence atoms as free particles in a vir- 
tually infinite potential environment to characterize their 
reactive behavior. In this regard, it is worth considering 
for such atoms the uniform motion by having dp/dt = 0 
through the time-constant associated wavefront condition 
and action S{r=cnst.,t)=cnst. (equivalent with Lagrangean 
constancy), in all given chemical space-points (atomic 
basins within molecule complex) [35]. This picture is also 
equivalently to have 



p(r) = if(r) = R 2 (r) 



(6) 



dS(r, t) 
dt 



0 



(12) 
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applied to eq. (8). By doing so, one obtains 



(VrS(r,*)r 
1m 



h 1 V 2 R{r) 
1m R(r) 



V(r) 



which can be rearranged as follows: 



-V n 



V(r) 



(13) 



(14) 



such that the total energy of a the valence system is now 
entirely driven by the quantum potential: 



V(r) 



■Vn 



(15) 



At this point, one can see that when turning to elec- 
tronegativity and combining eq. (15) with DFT definition 
(3), one obtains a generalization of the previous Boeyens 
formulation [6]: 



Xq-dft 



3Vq 
dN 



(16) 



V{r) 



which is the variation in the quantum potential with 
electron exchange under a constant classical or external 
potential. 

However, for a quantum characterization of the valence 
state, we are interested in how the energy described by 
eq. (15) varies under a quantum potential (10) 



Eq{o..u. 



V 2 R{r) V 2 p l l 2 {r 



1R{r) 1p l / 2 {r) 
lV 2 p(r) l[V rP (r)} 2 
4 



p(r) 



P 2 ( r ) 



(17) 



when the above relations (6) and (10) are substituted into 
eq. (15). 

It is worth noting that although we obtained the total 
energy (17) in the Bohmian mechanics context, it show- 
cases a clear electronic density dependency, not under a 
density functional (as DFT would require) but merely 
as a spatial function, which is a direct reflection of the 
entanglement behavior of Bohmian theory through the 
involvement of a quantum potential. However, in most 
cases, and especially for atomic systems, eq. (17) will 
yield numerical values under custom density function 
realizations. 

Golden ratio imbalance for valence states of atoms 

Atomic stability and periodicity remain major issues in 
the structural theories of matter; fortunately, they both 
have been largely solved by wave-particle (W7P) comple- 
mentarily quantum behavior; phenomenologically, such 
relationship can be expressed as "WAVE ® PARTICLE = 



constant", while it may be quantized (by Planck's con- 
stant h) in the light of Heisenberg principle as [36] 



WA VE&PAR TICLE = n w/P h 



(18a) 



Remarkably, when fixing the particle's observable 
property, say O, while letting wave information to vary, 
say AO, equation (18a) takes the workable form 



AO x 0 = n 0 h 



(18b) 



having as the preeminent realization the Bohr-de Broglie 
formulation a , leading with the first rationalization of the 
atomic periodicity [37]. However, when about the atomic 
chemical reactivity a similar analysis may be provided in 
terms of the number of electrons to atomic number 
ratio {N/Z): one may fix the observable ("particle") char- 
acter of the reactive atomic system by the ratio itself 



O 



N 
Z 



(19a) 



while modeling its evolving ("wave") character by the 
natural variation of the previous ratio in terms of 
exchanged electrons respecting the neutral state: 



AO 



AN N 



(19b) 



When combining eqs. (19a) and (19b) into eq. (18b) 
on the lowest quantized state (n 0 =l), the "ground state" 
of atomic reactivity that is the atom in its valence state 
so to speak, and within atomic units' formulation (i.e. by 
putting h=l, since the actual reactivity quantification 
involves only numbers with no dimension), one has the so 
called Heisenberg imbalance equation for valence atoms 



N-Z N 

x — = 1 

Z Z 



(20a) 



that can be rewritten as 

Z 2 + NZ - N 2 = 0 (20b) 
Eq. (20b) has the elementary acceptable solution 



—N + \/N 2 + 47V 2 

Z = = Nt 

1 



(21a) 



which establishes, the direct "chemical" connection be- 
tween the number of electrons and the atomic charge by 
means of the golden ratio 



-1 + V5 



0.6180 



(21b) 



generalizing the "physical" connection between nuclear 
(cosmic) synthesis at high pressure and atomic stability 



Putz Chemistry Central Journal 201 2, 6:1 35 
http://journal.chemistrycentral.eom/content/6/1/135 



Page 5 of 16 



in the gas phase (Z=A/); one has therefore the actual 
physical-to-chemical electronic charge - atomic number 
relationships 



1 . . .STABLE(PHYSICAL)A TOM 
N~ It... REACTIVE(CHEMICAL)ATOM 



(22) 



Worth remarking the results of type (20) and (22), here 
based on chemical reactivity specialization of Heisenberg 
type equations (18a) and/or (18b), were previously 
obtained at the level of neutron-protonic imbalance, in- 
side the atomic nuclei, based on well-founded empirical 
observations [6]. The present golden ratio appearance is 
ultimately sustained also by the deviation from the N=Z 
condition for so-called "quark atoms" (as another way 
in considering the atoms in a quantum valence state), 
earlier identified as true matter's entities responsible for 
matter's reactivity at the atomic level [38]. 

Therefore the atomic structure branching (22) can be 
regarded as the present golden ratio extension to valence 
atom and as such employed; actually, its consequences 
regarding the characterization of the quantum valence 
states of atoms within the Bohmian quantum potential 
are the main aims of the present endeavor and will be 
discussed next. 

Atomic implementation and discussion 

On Slater density for valence atoms 

Density is considered a "goldmine" in current computa- 
tional and conceptual quantum chemistry due to its link 
with observable quantities, energy density functionals in 
particular, as celebrated by DFT [13,20,39,40]. However, 
to quantitatively approach the chemical phenomenology 
presented in Figure 1, involving the ionization-to-affinity 
atomic description, the general Slater [33] density (in- 
volving the orbital parameter <f dependency) will be here 
employed for the first trial on modeling the combined 
Bohmian and gold-ratio features of valence atom; it 
assumes the general (trough still crude) working form: 



(23a) 



For the reactivity at the valence atomic level, or for 
some outer shell («) considered at the atomic frontier, 
one may assume almost electronic free motion or at 
least electronic motion under almost vanishing nuclear 
potential V(r); this way the density (23a), while entering 
the quantum potential (10) recovers the negative kinetic 
energy by the virial identity (14). Analytically, since eqs. 
(6), (10) and (23a), one has V 2 r p in = ?p V2 and the actual 
valence atomic virial realization looks like 



h 2 e 

2m 



2m 



(24) 



Equation (24) leaves with the identity: 

h£ = p 



(25a) 



that may be further rewritten with the help of the atomic 
Bohr-de Broglie relationship (see the note a ) to provide 
the atomic frontier radii shell-dependency 



^frontier 



(25b) 



Remarkably, the same result is obtained when employing 
a far more reach atomic shell structure description, namely 
when starting with the full atomic radial Schrodinger 
density [25] 



(2£) 



2k+1 



„2k-2 



(2«)! 



exp[— 2£ r] 



(26) 



and imposing the null-gradient condition [41], V r p n (r, f) = 
0, in accordance with the celebrated Bader condition of 
electronic flux of atoms-in-molecules [9,42], to yield: 



n 
1 



(25c) 



The identity between eqs. (25b) and (25c) gives suffi- 
cient support to the present Slater density approach 
eq. (23a) in modeling the valence atoms or the atoms 
at their frontiers approaching reactivity (i.e. atoms-in- 
molecules complexes by chemical reactions). 

Quantum chemical bonding and reactivity indices 

Once convinced by the usefulness of the Slater density 
form (23a) for the present valence atomic analysis, one 
will next employ it under the so called Parr-Bartolotti 
form [33] 



p(r,0=JV— exp[-2fr] 

7T 



(23b) 



such that to obey the ^-normalization condition, as 
required by DFT [43-47], 



47T r 2 p(r, Qdr = N 



(27) 



o 



by applying the Slater integral recipe 

oo 

f n\ 
I x" exp(—ax)dx = t 

o 



(28) 



It nevertheless showcases the parametric £ dependency 
that can be smeared out by considering the variational 
procedure 



0 



(29) 
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upon applying the total atomic energy 
£[fl = T[Z] + V m [Q + V„ e [Z) 



(30) 



where the components are individually evaluated within 
a radial atomic framework with the respective results for 
[21,48] 

• kinetic energy 



Anr 1 



i a 

r 2 dr V dr 



o 

e 

N — 
2 



nucleus -electronic interaction 



2r 



p(r, £)dr 



(31a) 



V„ e [fl = - j 471 r 2 = -AT£ (31b) 

o 

• inter-electronic interaction (see also Appendix) 



2N 



5 



(TV 2 -N)— £ 
v ; 16 



(31c) 



With these results, the optimum atomic parameter is 
quantified by the electronic number as follows: 



21 - 5N 
16 



(32) 



which immediately releases the working electronic 
density 



N (2\ - 5N 
n \ 16 



exp 



21 - 57V 



16 



2r 



(33) 



Having the completely analytical density in terms of 
number of reactive electrons as in eq. (33), worth point- 
ing here on the so called sign problem relating with its 
variation, e.g., its gradient, the gradient of its square 
root, etc.. Although this problem usually arises in density 
functional theory when specific energy functional are 
considered in gradient forms, see for instance ref. [49], 
there is quite instructive discussing the present behavior 
and its consequences. 

For instance, one can adapt either eqs. (25b) or (25c) 
through considering the present form (32) for the orbital 
exponent to be 



16 



n 21 - 57V 'bonding 



(25d) 



Here, one combines the frontier and maximum atomic 
radii with atoms-in-molecules phenomenology, as above 
indicated, to arrive to the present identification for the 
number of valence electrons possible to be involved in 
the same chemical bonding state as being Nbonding in 
(25d). Accordingly, the Figure 2 reveals interesting fea- 
tures of the present Slater-Parr-Bartolotti atomic density 
with quantum potential: 

• the fact that the (covalent) bond length is 
proportional to the atomic radii and in inverse 
correlation with bonding order is well known [50], 
and this it is also nicely reflected in eq. (25d); 
however, changing the sign to negative radii as 
surpassing the threshold 21/5 and fixing in fact the 
limit Nbonding=^> is consistent with maximum bond 
order met in Nature; it is also not surprising this 
self-released limit connects with golden ratio by the 
golden-spiral optimization of bond-order [51]; more 
subtle, it connects also with the 477 symmetry of two 
spherical valence atoms making a chemical bond 
(Figure 2, inset): such "spinning" reminds of the 
graviton symmetry [52] (the highest spherical 
symmetry in Nature, with spin equal 2) and justifies 
the recent treatments of chemical bonding by means 
of the quasi-particles known as bondons [10,53], as 
well as the use of the 4D complex projective 
geometry in modeling the chemical space as a non- 
Euclidian one, eventually with a time-space metrics 
including specific "gravitational effects" describing 
the bonding [51]; 

• the "gap" between the atomic systems contributing 
2 to 3 electrons to produce chemical bond is about 
double of the golden ratio, (r / '«) N bondi „ g =3 ~ 

i r / n ) N bondi „ g =2 — ^ T ' therefore, this gap marks the 
passage from the space occupied by a pair of 
electrons and that required when the third electron 
is added on the same bonding state: it means that 
the third electron practically needs one golden 
measure (r) to (covalently) share with each of the 
existing pairing electrons, while increasing the 
bond order to the level of three; it is therefore a 
space representation of the Pauli exclusion principle 
itself, an idea also earlier found in relation with 
dimensionless representation of a diatomic 
bonding energy (2r) at its equilibrium bonding 
distance r [54]; when the fourth electron is 
coming into the previous system, in order the 
maximum fourth order of bonding to be reach 
the chemical bonding space is inflating about 
five times more, yet forbidding further forced 
incoming electrons into the same space of 
bonding state as the bonding radius becomes 
negative in sign. 
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N, 



bo?iding 



Figure 2 Representation of the bonding length as a function of bonding electrons from valence atoms in molecule(s), based on eq. 
(25d), while marking the double golden ratio 2t gap between the bonding lengths of the second and third bonding order, as well as 
the forbidden chemical bonding region for N bonding >21/5 for the electrons participating in the same bonding. Further connection of 
chemical bonding and the 4D space to model it is suggested by the inset picture illustrating the 2-fold (477) spinning symmetry of the adduct 
atom respecting the bonding direction, after [51]. 



Having revealed the chemical bonding information car- 
ried by the density (33) when considered for combined 
valence atoms-in-molecules, it is next employed on ener- 
getically describing the atomic reactivity as a propensity 
for allowing electronic exchanging and bonding. As such, 
it leaves the total quantum (Bohmian) energy in (17) with 
the compact form 



Eq\{o..U. 



25 . 

AT 

512 

5^(21 -5A0 2 



105 441 

N + 

256 512 



(34a) 



Note that the actual working total energy is not that 
obtained by replacing the density (33) in eqs. (31a)-(31c) 
and then in total energy (30) because here the double- 
variational procedure was considered; that is, the first 
optimization condition was considered as in eq. (29), and 
the resulting (optimum) density (33) was then employed 
in the quantum energy (17), which in turn was obtained 
by applying the variational eq. (12) to the perceived phase 
transition in the Bohm eikonal wave-function (5). To 
emphasize the accuracy of eq. (17) over that of (30) with 
density (33), when one considers the last case, eq. (30) 
yields the following non-quadratic form for energy: 



E\{a.u. 



25 

512 
N 

512 



-N 



;N 



105 . 
256 J 
(21 - 5N) 2 



441 
512 J 



-N 



(34b) 



which is not appropriate for describing the valence state 
of an atom, as eq. (2) prescribes, despite being similar in 
form to the Bohmian-based result of eq. (34a). Thus, the 
previous limitation of the Parr-Bartolotti conclusion [33] 
and the paradox raised in describing the valence (para- 
bolically) state with the optimized atomic density (33) are 
here solved by the double (or the orthogonal) variational 
implementation, as recently proved to be customary for 
chemical spaces [32] . In the light of this remark one may 
explain also the sign difference between the "physical" 
energy (34b) and that obtained for the "chemical" situ- 
ation (34a): through simple variational procedure for 
"physical" energy (30) the result (34b) is inherently nega- 
tive - modeling systems stability in agreement with the 
upper branch of eq. (22), whereas the double variational 
algorithm employing optimized density (33) into the 
Bohmian shaped energy (17) it produces the positive out- 
put (34a) associated with activation energy characteristic 
for chemical reactivity corresponding to the lower branch 
ofeq. (22). 

Therefore, to be accurate, one should consider the 
quantum potential related optimized energy (34a) instead 
of simply the orbital optimized one of eq. (34b). There- 
fore, assuming that eq. (34a) appropriately describes the 
atomic valence state in DFT (see the upper/reactive 
curve in Figure lb), the next task is to search for the 
quantum valence charge for which the valence energy 
approaches its optimum value (or the "ground state" of 
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the atomic chemical-reactivity, i.e. the previously golden- 
ratio quantification of the valence atomic state); to this 
aim, at this point, one can employ the golden ratio rela- 
tionship (21a) and first rewrite eq. (34a) as 

£qi(«.m.) =— (6.80r-A0 2 (35) 
oil 

which is minimized at the value 

N = 6.80r (36) 

However, one must again apply the double-variational 
procedure, now in terms of number of electrons, i.e., 
reconsidering eq. (36) with the golden ratio at the react- 
ive (chemical) electronic level of eq. (22) such that a sec- 
ond equation is formed 

JV= 6.80- (37) 
N 

with the positive solution 

N REACT =N** = 2.60768VZ" (38) 

This expression avails of the significance of the max- 
imum number of electrons, for a given atom, possibly 
engaged in a reactive environment by either (or both) 
accepting or (and) ceding electrons to or from its 
valence state, see Table 1. 

The result of this process is different from the expected 
physical result {N STABIL =Z) according to the upper 
branch of eq. (22), which is higher than the physical one 
until reaching the carbon system {Z INTRCHANGE =6.8), 
while continuing below it thereafter (see Figure 3). 

The above interchange (effective) atomic number 
through which the chemical (reactive) state is associated 
with lower charge respecting the physical state may be 
also be found at the energetic level based on quantum 
equation (34a), as specialized for the two branches of 
Figure 3 for the N(Z) dependence. Thus, the chemical 
(reactive) state takes the analytic form 

Eqi (n react — 2.60768\/z) 

= 0.861328 - 1. 06956 VZ + 0.332031Z (39) 

and interchanges with the ground state E q1 (N S table — ► 2) 
at the points {3.5,6.8}, as observed also from Figure 4; 
however, the interchanging point beyond which all 
chemical atomic systems are more stable in the chem- 
ical or reactive state than in the physical ground state 
is consistently recovered. 



Table 1 Synopsis of the critical charges in the physical 
ground state (N ) as well as for chemical reactive 
(valence) state (N**) for atoms of the first four periods of 
the periodic table of elements, as computed from the 



minimum point of associated interpolations of ionization 
and electronic affinities [33] and of eq. (38), respectively 



Atom 


Z 


l[ev] 


A[ev] 


N* 


N** 


H 


1 


13.595 


0.7542 


0.558735 


0.607681 


Li 


3 


5.390 


0.620 


0.629979 


0.516636 


B 


5 


8.296 


0.278 


0.534672 


0.830952 


C 


6 


11.256 


1.268 


0.626952 


0.387488 


0 


8 


13.614 


1.462 


0.620309 


0.375636 


F 


9 


17.42 


3.399 


0.742422 


0.823043 


Na 


11 


5.138 


0.546 


0.618902 


0.648699 


Al 


13 


5.984 


0.442 


0.579755 


0.402127 


Si 


14 


8.149 


1.385 


0.70476 


0.757049 


P 


15 


1 0.484 


0.7464 


0.576651 


0.0995049 


S 


16 


10.357 


2.0772 


0.750876 


0.430724 


CI 


17 


13.01 


3.615 


0.884779 


0.751744 


K 


19 


4.339 


0.5012 


0.630596 


0.366618 


V 


23 


6.74 


0.526 


0.584648 


0.505999 


Cr 


24 


6.763 


0.667 


0.609416 


0.774976 


Fe 


26 


7.90 


0.164 


0.5212 


0.296616 


Co 


27 


7.86 


0.662 


0.59197 


0.549908 


Ni 


28 


7.633 


1.157 


0.67866 


0.798551 


Cu 


29 


7.724 


1.226 


0.688673 


0.0427917 


Se 


34 


9.75 


2.0206 


0.761417 


0.205262 


Br 


35 


11.84 


3.364 


0.896885 


0.427249 


Rb 


37 


4.176 


0.4860 


0.631707 


0.861 904 


Zr 


40 


6.84 


0.427 


0.566584 


0.492423 


Nb 


41 


6.88 


0.894 


0.649348 


0.697305 


Mo 


42 


7.10 


0.747 


0.617582 


0.899704 


Rh 


45 


7.46 


1.138 


0.680006 


0.492856 


Pd 


46 


8.33 


0.558 


0.571796 


0.686153 


Ag 


47 


7.574 


1.303 


0.707782 


0.87736 


Sn 


50 


7.342 


1.25 


0.705187 


0.439089 


Sb 


51 


8.639 


1.05 


0.638358 


0.622567 


Te 


52 


9.01 


1 .9708 


0.779975 


0.804255 


1 


53 


1 0.454 


3.061 


0.91404 


0.984204 



Nevertheless, the energetic analysis also reveals the 
atomic systems Be, B and C to be situated over the cor- 
responding physical stable states; this may explain why 
boron and carbon present special chemical phenomen- 
ology (e.g., triple electronic bonds and nanosystems with 
long C-bindings, respectively), which is not entirely 
explained by ordinary physical atomic paradigms [55-60]. 

The energetic discourse may be complete with the 
electronegativity and chemical hardness evaluations by 
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Nreact -2.61VZ 




^'stable ~Z 




^■'stable ~ Nreact 




Figure 3 The comparative shapes of the valence electrons to 
be engaged in chemical reactivity (continuous curve) 
computed using eq. (38) based on the combined optimal Bohm 
total energy (35) with the golden ratio imbalance of eq. (22), 
respecting the stable physical case (dot-dashed curve), and of 
their differences (dashed curve); all originate at the 0 th atom 
(the neutron, Z=0). 



applying the DFT definitions (3) and (4) to physical and 
chemical energies, respectively. In the first case, expres- 
sion (34b) is applied to provide the following so-called 
"universal" forms of Parr and Bartolotti [33]: 



XpBl 



1PB1 



3£i 
dN 

1 fd 2 E l 



2\dN 2 



16 

135 



512 



(a.u.) = 5.l(eV) (40) 
{a.u.) = 7.17(eV) (41) 



The result, nevertheless, appears to be an unusually 
higher increase in chemical hardness than in electro- 
negativity, which certainly cannot be used to model a 
reactive-engaged tendency because it is more stable (by 
chemical hardness) than reactive (by electronegativity); it 
is, however, consistent with the physical stability of the 
system, provided by the single variational procedure 
through which eq. (34b) was produced. 

Instead, to chemically model reactivity, the double 
variation procedure is applied and eq. (34a) is substi- 
tuted into eqs. (3) and (4), though by considering also 
the double reactive procedure for charge as well, i.e., by 
considering eq. (38) with the golden ratio information of 
(22) to respectively yield the results 



Xpbi 




Af=2.60768v^/ z=rA7 / N=1 

0.209963(a.w.) = 5.713(eV) 



(42) 



1PB2 



a 2 £ 



Qi 



dN 2 



25 
512 



{a.u.) = 1.3286(eV) (43) 



Remarkably, the actual electronegativity of (42) obtained 
by the quantum Bohm and golden ratio double pro- 
cedure yields sensible results similar to those of the 
single variational approach (40); however, the chemical 
hardness of (43) is approximately 5-fold lower than its 
"stable" counterpart (41), affirming therefore the mani- 
festly reactive framework it produces - one described by 
a quadratic equation (34a) instead of a cubic one (34b). 

Charge waves in gauge chemical reactivity 

Finally, one considers the chemical reactivity discussion 
as based on the gauge reaction that equilibrates the 
chemical bond by symmetrical bond polarities [25] 



B + = A — B = A + +B~ 



(44a) 



such that the reactive electrons are varied on the 
reunited intervals of eq. (1); such analysis was previously 
employed to fundament systematic electronegativity and 
chemical hardness definitions by the averaging (through 
the integration) factor 



0.5 



1 



(44b) 



dN 



along the reaction path accounting for the acidic (electron 
accepting, 0 < N < +1) and basic (electron donating, 
-1 < N < 0) chemical behaviors. 

In this scaled (gauge) context of reactivity, the fore- 
going discussion is dedicated to investigating the link be- 
tween the critical ground state charge (N ) and the valence 
or reactive state (N ). While the first appears as a conse- 
quence of naturally fitting the three points in Figure 1 (the 
ionization, neutral and affinity states), with the effect of 
biasing the minimum of the energetic curve in Figure lb 
with respect to the apparent Parr-DFT curve in Figure la, 
and is thus derived graphically (see Figure 5), the valence 
charge is based on the combined quantum energy and 
golden ratio information in eq. (38). Both are reported 
for the indicated number of atomic systems of the peri- 
odic table of elements in Table 1. One notes, for instance, 
that while the critical ground state charge N always lies 
in the range [0.5,1], the valence charge N may span the 
interval [0,1]; one may interpret such behavior as being 
associated with the difference between the fraction % and 
integer "1" in driving the principles of chemical reactivity 
and the electrophilicity equalization principle in particu- 
lar, when the "quantum transition" 1/2— >1 is required in 
the energy exchange of chemical systems for it to be valid 
for both electronegativity and chemical reactivity prin- 
ciples [61]; nevertheless such scaling it is equivalent with 
above acidic-basic gauge averaging of eq. (44b). This 
way, the valence charge problem may be extended to the 



Putz Chemistry Central Journal 201 2, 6:1 35 
http://journal.chemistrycentral.eom/content/6/1/135 



Page 10 of 16 



E[a.u.] 






Eqreact 
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Figure 4 The same comparative shapes shown in Figure 3, here at the level of energy (34a) specialized for the reactive and the stable 
N(Z) dependencies of Figure 3; the various plots successively display increasingly large atomic Z-ranges to better emphasize the 
chemical vs. physical behavior (see text). 
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Figure 5 A graphical interpolation for selected elements of Table 1 in terms of their ionization, neutral and affinity states, aiming to 
determine the critical (displaced) charge of the DFT ground state, as prescribed by Figure lb. 



interval [0,2], at its turn seen as a gauge transformation 
of the chemical reactivity charge domain f-1, +1], where 
one reencounters the challenging problem of whether 
the "One electron is less than half what an electron pair 
is" [62], the response to which is generally complex but 
may here be approached through the following steps. 

First, by employing the data presented in Table 1, one 
constructs the so-called "continuous" ground and 
valence charge states by appropriately fitting over the 
first four periods of elements, here restrained to 10 th - 
order polynomials. This is performed by interpolating 
every three points of the 32 elements presented in 
Table 1, although by spanning the atomic number range 
Ze [1,53], thus yielding (see also the allied representa- 
tions of Figure 6): 



NC = 0.677771 - 0.193006Z + 0.104303Z 2 

- 0.0242757Z 3 + 0.00302359Z 4 

- 0.000220204Z 5 + 9.79371- lO^Z 6 

- 2.6894-10" 7 Z 7 + 4.4448-10" 9 Z 8 



4.05002 10 _li Z y 



1.56254- 10~ 13 Z 10 (45a) 



NC" = 0.768074 - 0.224502Z + 0.076654Z 2 

- 0.0107337Z 3 + 0.000667575Z 4 

- 0.000012746Z 5 - 6.0246 10"^ 



4.01411 10" 8 Z 7 



9.49568- 10" 10 Z 8 



1.05449 10" n Z 9 - 4.58533 10" 14 Z 10 



(45b) 



Putz Chemistry Central Journal 201 2, 6:1 35 
http://journal.chemistrycentral.eom/content/6/1/135 



Page 12 of 16 



N" 




X 

1.2 
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Figure 6 The critical ground state and valence charge points for the elements of Table 1 and their 10 th -order continuous 
interpolations according to eqs. (45a) and (45b). 



Equations (45a) and (45b) are then combined into a The result (47) has the following conceptual fundamen- 
sort of special charge wave function based on their dif- tal quantitative interpretation: the difference between the 
ference on the golden ratio scale (see Figure 6 for graph- ground and valence optimum charges is regulated by the 



ical representation) 

W z = t(NC* z - NC* Z 



(46) 



with the peculiar property that its square-integrated 
form over the Z-range of interpolation gives 



I W z \ dZ= 0.667233 = r 



golden ratio scale, or in other terms, 



53 



I (NC* Z - NC* z *fdZ « - = 1 + t 



(48) 



such that it provides a sort of normalization corrected by 
(47) the golden ratio value; it also fulfills the interesting rela- 
tionship: 
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\ 



53 



(JVC* 



NC*z) dZ<=ar 



(49) 



In any case, the present analysis provides the qualita- 
tive result that the difference between the critical ground 
state and optimal valence charges is more than half of 
an electronic pair, giving rise to the significant notion 
that chemical reactivity is not necessarily governed by a 
pair of electrons but governed by no less than half of a 
pair and is related to the golden ratio (r > 0.5). 

However, fractional values in general and those related 
to the golden ratio particular, may be interpreted as a 



consistent manifestation of the quantum mechanical (i.e., 
wave functional) approach of chemical phenomena, here 
at the reactivity level. Moreover, the quadratic critical 
charge function (46), as shown in Figure 7, clearly reveals 
that a higher contribution to electronic pair chemistry 
is given by the third period of elements and by the 
third and fourth transitional elements in particular, a 
result that nicely agrees with the geometrical interpret- 
ation of the chemical bond, particularly the crystal lig- 
and field paradigm of inorganic chemistry [9]. 

Also a local analysis of the type of charge that is 
dominant in atomic stability, i.e., the critical physical 
ground state or the chemical valence reactive state 
based on eqs. (45a) and (45b), respectively, may be of 
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considerable utility in refined inorganic chemistry 
structure-reactivity analysis. To the same extent, it 
depends on the degree of the polynomials used to 
interpolate the critical and valence charges over the con- 
cerned systems; however, through the present endeavor, 
we may assert that the analysis should be of the type 
(48), which in turn remains a sort of integral version of 
the imbalance equation (20a), in this case for the 
ground-valence charge gap states of a chemical system. 

Conclusions 

Aiming to hint at the solution to the current debate 
regarding the physical vs. chemical definition of an 
atom and as a special stage of a larger project regard- 
ing quantum chemical orthogonal spaces, the present 
work addresses the challenging problem of defining 
and characterizing valence states with respect to the 
ground state within conceptual density functional the- 
ory. We are aware of the earlier warnings raised by Parr 
and Bartolotti and others [18,33,63] regarding the limits 
of density functional theory and of the total energy of 
atomic systems combined with a Slater-based working 
density to provide a quadratic form in terms of system 
charge, as required by the general theory of chemical re- 
activity of atoms and molecules in terms of electronega- 
tivity and chemical hardness. Fortunately, we discovered 
that the Bohmian form of the total energy of such atomic 
systems provides, instead, the correct behavior, although 
it is only density-function-dependent and not a func- 
tional expression. Moreover, this finding was reached 
through the so-called double variational procedure, 
which, as emphasized earlier, was likely to reproduce 
the chemical reactivity principles of electronegativity and 
chemical hardness in an analytical manner; however, 
such a double analytical variational approach is consist- 
ent with the recent advanced chemical orthogonal spaces 
approaches of chemical phenomenology [64] as being at 
least complementary to the physical description of many- 
electronic systems when they are engaging in reactivity 
or equilibrium as the atoms-in-molecules Bader theory 
prescribes [9,42]. With the present Bohmian approach, 
the total energy is in fact identified with the quantum 
potential, thus inherently possessing non-locality and 
appropriate reactivity features, which are manifested 
even over long distances [10,11,53]; this also generalizes 
the previous Boeyens electronegativity formulation of 
electronegativity [5,6] from the direct relationship be- 
tween a quantum potential and its charge derivative. The 
double algorithm was also implemented to discriminate 
the valence from the ground state charges, this time by 
using the golden ratio imbalance equation as provided by 
adaption of the Heisenberg type relationship to chemical 
reactivity for atoms. This corresponds to an analytical 



unfolding of the physical and chemical imbalance of the 
electronic charge stability of atomic systems, paralleling 
the deviation from the equal electron-to-proton occu- 
pancy in physical systems toward electron deficiency in 
the valence states of chemical systems. This dichotomy 
was implemented by the golden ratio presented in eq. 
(22). As a consequence, the difference between valence 
and ground state charge systems is naturally revealed 
and allows for the explanation of chemical reactivity and 
bonding in terms of fractional electron pairs, althrough 
driven by the golden ratio under the so-called physical- 
to-chemical charge difference wave function and asso- 
ciated normalizations, all of which represent elaborated 
or integral forms of the basic imbalance atomic equation. 
The present results are based on 10 th -order polynomial 
fitted over 32 elements from the first 54 elements of the 
first four periods of periodic table of elements and can be 
further pursued by performing such systematic interpola- 
tions that preserve the golden ratio relationships, as 
advanced herein; they may also provide a comprehensive 
picture of how valence electrons may always be pro- 
jected/equalized/transposed into ground state electrons 
within the perspective of further modeling chemical reac- 
tions when chemical reactivity negotiates the physical 
molecular stabilization of atoms in molecules. 

Endnotes 

a For circular orbits, the lowest ones in each atomic 
shells - including the valence ones, one has AO=Ar=2nr, 
with r the orbital radii thereof, while 0=p is the fixed 
particle's momentum on that orbit; therefore, when 
combined into eq. (18b) they provide the celebrated 
Bohr-de Broglie relationship rp=nh solving the atomic 
spectra of Hydrogen atom in principal quantum num- 
bers («). 

Appendix: Semi-classical inter-electronic energy 

For the inter-electronic interaction, see Figure 8; in 
evaluating V ee [<f] of eq. (31c), the two-electronic density 
is approximated by the Coulombic two mono-electronic 
density product, thus neglecting the second-order density 
matrix effects associated with the exchange-correlation 
density. 

However, for the analytical evaluation of the electron- 
electron repulsion energy using the density (23b), much 
care must be taken. For instance, one has to use the 
electrostatic Gauss theorem, which states that the clas- 
sical electrostatic potential outside a uniform spherical 
shell of charge is just what it would be if that charge 
were localized at the center of the shell and that the po- 
tential everywhere inside such a shell is that at the sur- 
face, [21,48] see Figure 8. Therefore, the electronic 
repulsion energy becomes 
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which recovers the expression presented by eq. (31c), 
when the Slater integral type of Eq. (28) is also employed. 
Note that the electron-electron repulsion term was 




Figure 8 Representation of the space regions of the 1 and 2" 
electrons, their potential influences and reciprocal interaction 
[21,48]. 



written by also considering the Fermi- Amaldi (N-l)/N 
factor [13], which ensures the correct self-interaction be- 
havior: when only one electron is considered, the self- 
interaction energy must be zero, V ee (AT— >1)— >0. 
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